Improved efficiency with variational Monte Carlo using two level sampling 
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A two level sampling method is applied to variational Monte Carlo (VMC) that samples the one 
and two body parts of the wave function separately. The method is demonstrated on a single Li2 
molecule in free space and 32 H2 molecules in a periodic box. This simple modification increases 
the efhciency of a VMC simulation by up to 72%. 



o 
o 

O ■ 

Oh: 



> 
in 

o 
l> 
o 
o 
o 

o 

• 1— I 
Oh 



X 



Variational Monte Carlo (VMC) § is an effective 
method for computing the ground state properties of 
atoms, |, I molecules, |, | and solids. [| 0, |, | Ex- 
plicitly including the electron correlation into the wave 
function allows VMC to recover a large percentage of the 
correlation energy, typically 70-90% (or more for small 
atoms and molecules) . ||, ^, [To| 

VMC computes matrix elements with respect to the 
square of a trial wave function, ^p. The most common 
matrix element is the energy. 



E = 



J dR ^2(i?; a)EL{R; a) 



(1) 



where El = ^Hip. The variational principle states that 
this energy will be greater than or equal to the true 
ground state energy for any tJj. Typically, the wave func- 
tion is parameterized and those parameters are varied 
until a minimum energy |l0[ (or alternatively, minimum 
variance [§, ^) is reached. 

The Metropolis method ||ll|] constructs a transition 
probability, P(s' — > s), for use in a Markov process. The 
result of that process is the desired normalized probabil- 
ity distribution. A correct transition probability can be 
constructed by satisfying detailed balance, 



Tr{s)V{s s') = Tr{s')V{s' s) 



(2) 



where 7r(s) is the desired probability distribution {tp'^{R) 
for VMC ). In the generalized Metropolis method, the 
transition probability is split into two pieces, an a pri- 
ori sampling distribution, T(s — > s') and an acceptance 
probability A{s s'). 

The original Metropolis method generates a trial move, 
R' , of one particle inside a box of side A centered about 
the old position, R. This leads to a constant T inside the 
box and zero outside. The trial move is then accepted or 
rejected with probability 



A 



1, 



iP^iR) 



(3) 



In electronic problems, a typical wave function is 

iP = De-'^ (4) 

where D is the product of spin up and spin down Slater 
determinants of single particle orbitals and U is the two 
(or higher) body Jastrow factor. This latter part contains 
the explicit electron-electron correlation. 

A simulation step consists of a trial move of every elec- 
tron. Each trial move is accepted or rejected according 
toEqd). 

A single electron is moved at a time, which only 
changes one column of the Slater matrix. By saving and 
using the inverse of this matrix, computing the new de- 
terminant and its derivatives is an 0{N) operation. 
If the move is accepted, the inverse matrix can be up- 
dated with 0(7V^) operations. (Note that computing the 
determinant directly is of order 0{N^).) A smaller ac- 
ceptance ratio will be faster, since fewer updates need to 
be performed. 

Multilevel sampling has been used extensively in path 
integral Monte Carlo. ||l^ In multilevel sampling, the 
wave function is split into several factors (levels), and 
an accept /reject step is performed after computing each 
factor. The entire wave function is computed only if there 
is an acceptance after every factor. Rejections will make 
the overall algorithm faster, since not all the factors need 
to be computed. For VMC it is natural to use two factors 
- the single body part, D, as one factor and the two body 
part, e~'^ , as the other. The single body part is the 
cheaper of the two to compute, so computing it first will 
give the greatest reduction in time. 

The two level sampling algorithm for VMC proceeds 
as follows. A trial move, R' , is proposed and accepted at 
the first stage with probability 



D^R') 



D^R) 



(5) 



If accepted at the first stage, the two body part is 
computed and the trial move is accepted with probability 



A2 = 



1, 



exp [-2f7(i?')] 



exp [-2C/(i?)] 



(6) 
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It can be verified by substitution that these satisfy de- 
tailed balance in Eq. (0). After an acceptance at this 
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second level, the inverse Slater matrices are updated as 
described previously. 

The figure of merit for these simulations is the effi- 
ciency 



1 



(7) 



where T is the computer time and a is the error estimate 
of the average of some property such as the energy. 

Maximum efficiency results from a competition be- 
tween two effects. The first is that a lower acceptance 
ratio will yield a lower run time and a larger efficiency. 
The second is serial correlations. A lower acceptance ra- 
tio will give a longer correlation time, hence a larger error 
and lower efficiency. 

In maximizing the efficiency, the obvious parameter to 
adjust is the step size, A. But the number of steps be- 
tween computations of El can also be adjusted. The 
Metropolis algorithm produces correlated state points, 
so successive samples of E]^ don't contain much new in- 
formation. It is advantageous to compute E^ every few 
steps rather than every step. In this work the local en- 
ergy was computed every 5 steps. 

The algorithm is demonstrated using a Li2 molecule in 
free space and a collection of 32 H2 molecules in a peri- 



odic box. The wave functions, which are the ^'m's from 
Reynolds, et ai, use a simple electron-electron and 
electron-nuclear Jastrow term, and use floating Gaus- 
sians for the orbitals. The hydrogen molecules were in 
box of side 19.344 atomic units (r^ = 3.0). 

The results for the different sampling methods for the 
Li2 molecule are given in Tables | and |n[ The second 
level acceptance ratio is quite high, indicating the single 
body part is a good approximation to the whole wave 
function. The efficiency is also shown in Figure 1. 

The results for the different sampling methods for the 
H2 molecules are given in Tables HI and IV. The effi- 
ciency is also shown in Figure 2. 

Comparing the maximum efficiency for each sampling 
method, two level sampling is 39% more efficient than 
the standard sampling for Li2, and 72% more efficient 
for 32 Hz's. 

More complicated schemes using more levels or differ- 
ent splittings of the wave function could be devised. This 
particular scheme is attractive because it uses quantities 
readily available in a VMC computation, and requires 
minimal modification to existing VMC sampling algo- 
rithms. 

This work has been supported by the computational 
facilities at NCSA and by NSF grant DMR 98-02373. 
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TABLE I: Timings for Li2 molecule using the standard sam- 
pling method. All times in seconds on an SGI Origin 2000. 
Acceptance Determinant Jastrow Total 



A 


Ratio 


Time 


Time 


Time 




1.0 


0.610 


48.3 


340 


516 


1190 


1.5 


0.491 


48.1 


340 


508 


1680 


2.0 


0.407 


48.2 


340 


503 


1460 


2.5 


0.349 


48.2 


339 


499 


1070 


3.0 


0.307 


48.2 


339 


496 


800 



TABLE II: Timings for Li2 molecule using the two level sam- 
pling method. All times in seconds on an SGI Origin 2000. 



A 


First Level Second Level Total Acc. 
Acc. Ratio Acc. Ratio Ratio 


Time 




1.0 


0.674 


0.899 


0.606 


400 


1580 


1.5 


0.543 


0.894 


0.485 


347 


2430 


2.0 


0.447 


0.897 


0.401 


304 


2340 


2.5 


0.379 


0.902 


0.342 


276 


1910 


3.0 


0.331 


0.906 


0.300 


256 


1400 



TABLE III: Timings for 32 H2 molecules in a periodic box 
using the standard sampling method. All times in seconds on 
a Sun Ultra 5. 



Acceptance Determinant Jastrow Total 



A 


Ratio 


Time 


Time 


Time 




2.0 


0.606 


167 


1089 


2015 


0.61 


3.0 


0.455 


167 


1085 


1891 


1.22 


4.0 


0.338 


166 


1084 


1794 


1.23 


5.0 


0.250 


166 


1080 


1722 


1.06 


6.0 


0.185 


164 


1080 


1668 


1.02 


7.0 


0.139 


162 


1084 


1629 


0.76 



TABLE IV: Timings for 32 H2 molecules in a periodic box 
using the two level sampling method. All times in seconds on 
a Sun Ultra 5. 

First Level Second Level Total Acc. Total 



A 


Acc. Ratio 


Acc. Ratio 


Ratio 


Time 




2.0 


0.740 


0.795 


0.589 


1804 


0.59 


3.0 


0.598 


0.728 


0.436 


1421 


1.77 


4.0 


0.468 


0.681 


0.319 


1185 


2.11 


5.0 


0.357 


0.649 


0.232 


994 


1.55 


6.0 


0.370 


0.627 


0.169 


849 


1.87 


7.0 


0.204 


0.609 


0.124 


740 


1.46 





FIG. 2: Efficiency of VMC for 32 H2 molecules. 



